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Abstract 
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1 Introduction 



In this paper we improve constraints on the masses of the scalar excitations of the minimal 
super symmetric Standard Model (MSSM) (for a review of the MSSM see ref. |]T|) by 
making use of a new method of computing vacuum tunneling rates in field theories 

To begin with, recall that the most general potential for the scalar fields in the MSSM 
that respects the gauge symetries and conserves baryon and lepton numbers is: 

v = E 

3 

+ J2 m %^3 + AuhuQ L H 2 U R + A D h D Q L E x b R 

hi 

+ AeJieLlHiEr + BiiHiH 2 

where W is the superpotential 

W = huQ L H 2 U R + h D Q L H x D R + h E L L H x E R + [ iE x H 2 (1.2) 

and 4>j stands for any scalar field in the theory (rrijj are gauge invariant mass terms). 
Hi and H 2 are the Higgs superfields, while Ql, U r , and D R stand for up and down type 
quark superfields. L L and E R are the lepton superfields (Ql and L L are weak SU(2) 
doublets). The Yukawa couplings h v , h D , h E , as well as the products Ajjhu, A D h D , and 
A e Ue are 3x3 matrices in generation space. Superpartners are denoted by a tilde over the 
corresponding superfield. The last five terms in equation ( |1 . 1| ) are the soft supersymmetry 
breaking terms. The expression for the D term is 

V D = \E9l{E& a ^) (1-3) 

where g a (a = 1,2,3) are the gauge coupling constants and T a are the generators of the 
gauge group algebra in the representation of the scalar fields (p a (namely, the squark, 
slepton and Higgs fields). 

Even at the classical level, the potential in equation ( |1 . 1| ) is a complicated quartic 
function which has many local minima. A systematic study reveals that for specific 
ranges of parameters the potential is unbounded from below (UFB). Several charge and/or 
color breaking (CCB) minima of the potential may also be present || |], ||. Later we will 
recollect some of the constraints on the parameter space arising from the requirement that 
the charge and color preserving minimum be a global minimum of the potential. This 
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requirement is however too conservative. Indeed, it is quite possible that the vacuum 
we live in is a metastable false vacuum with a very long life-time. There have been 
some numerical studies of vacuum decay rates in restricted sectors of the MSSM |6|, 
f7|. However, untill recently no techniques existed for efficiently scanning the parameter 
space for unstable vacua in the full theory because of the large number of fields involved. 
Recently a new method for simplifying the computation of vacuum tunneling rates in field 
theories with many scalar fields was introduced by one of us @ (in this context, see also 
ref. f8|). The goal of the present work is to use this method to relax the mass constraints. 
The constraints we will obtain by computing a lower bound on the tunneling amplitude of 
the viable MSSM vacuum will be "necessary conditions" and the parameter space ruled 
out by these constraints is "necessarily" ruled out. 

Our numerical method is described in section 2. In section 3 we use this method to 
constrain the allowed ranges of particle masses. 



2 The Method of Reduction 



Suppose we wish to find vacuum tunneling amplitudes in a 4- dimensional field theory 
described by the (Euclidean) action 



S 
T 



T + V 

d 4 x 



(2.4) 



2 r 



dx. 



V = J d 4 x[U(<P U ..4n)] 



where (f)j, j = 1, 2..., n are real scalar fields and U is a potential. Let us also suppose that 
U has a local minimum at (the false vacuum) and a global minimum at </>*■ (the true 
vacuum). We normalize U so that U{m) = 0. Then the probability density for the rate 
of nucleation of true-vacuum bubbles in the false vacuum is given by 



r 

V 



Aexp 



h 



[i + 0(h)] 



(2.5) 



where <pj is a saddle point of the action S (known as the bounce) and A is a pre-factor of 



dimension (mass) 4 || nj [11]]. Let us summarize salient features of equation (|2.5|) : 

(i) The bounce is a finite action field configuration which is invariant under time-reversal. 



It follows that (j>j(x) 



<pt as x 



oo and d(pj(x)/dt \ t - = 0. 
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(ii) Equation (|2.5|) is the field theoretic version of a generalized semiclassical WKB tun- 



neling formula in quantum mechanics fll2| . The field theoretic generalization involves 
renormalization of both S and A ||11|| . If we choose a renormalization scheme where the 
masses and couplings entering in S are the physical masses and runing couplings at some 
renormalization scale Q } then and A are both functions of Q in a way such that y is 
Q independent. For the MSSM, it has been argued that as long as the mass parameters 
and the inverse size of the bounce <pj are within one or two orders of magnitude from 
the electroweak scale Qo, the value of — j^- and therefore A is relatively insensitive to 
radiative corrections provided Q is chosen to be of order of the electroweak scale J|, [7[]. 
In practice, we will choose Q = Qq ~ 200GeV, find the bounce, compute S and make 



a simple estimation of A which is calculable in principle |]TTj, |13| but is difficult to cal- 
culate exactly even in very simple field theories. We will write A as rjv 4 where v is the 
characteristic mass scale of the theory (in the present case v ~ Qo ~ 200 GeV) and 77 is 
a factor representing our ignorance. Smallness of the radiative corrections imply that the 
error involved in making the approximation rj ~ 1 is negligible (except for a possible en- 
hancement described below). The approximations made above are valid when the scalar 
masses and the inverse of bounce size are not off by several orders of magnitude from Q . 

(iii) The probability of vacuum decay is large if T/V x TJ 4 ~ 1 where To is the age of 
the universe. The condition T/V > Tj 4 translates roughly into 

^Sl < 400 (2.6) 

when 77 is taken to be of order unity. Note that close to the surface separating the regions 
in the parameter space where (T/V x Tq) — 1 changes sign, a 5% decrease in the tunneling 
rate would require 77 to be as large as exp(20). There is a recent claim that 77 may be 
enhanced by a factor of [^IM]^/ 2 if the bounce breaks a continuous internal symmetry 
group of the false vacuum with N generators |14[]. We however choose to be conservative 
and do not apply this correction which can only strengthen the bounds we obtain on 
the parameter space in the end. Also note that we consider only a zero temparature 
quantum tunneling effect. There is a corresponding finite temperature effect for which 
the tunneling rate will be larger but the available time for transition will be smaller. The 



present analysis may be extended along the lines of [0 to obtain independent constraints 
from the finite temperature case. 

The most important part of the calculation is the determination of S[<j)j] for which we 
use the following method of reduction (more details can be found in ref. ||). Recall that 
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the bounce is a saddle point of the action with a single unstable direction. It satisfies 
stationarity with respect to scale transformations x — *> Ax, where A a positive number. In 
a d dimensional field theory (d > 2) the condition of stationarity is |16| 



(d-2)T+(d)V = 0. (2.7) 

An interesting property of the bounce is that it is the global minimum of the action in 
the subspace C of all field configurations with the trivial boundary condition that also 
satisfy equation (|2.7|). A rigorous upper bound on the bounce action is obtained by 
taking any field configuration (j>j(x) having trivial boundary conditions for which V < 
and scale transforming it to satisfy equation ( p.7|) . A better upper bound is obtained by 
minimizing the action on a finite dimensional subspace of C. A good numerical algorithm 
should start with a suitably chosen finite dimensional space over which the minimization 
is a quick procedure and pays off well in lowering the action. This is the objective of our 
computations, namely, to obtain the best possible lower bound on the tunneling rate (given 
time and computer limitations) and use it to rule out as much of the parameter space as 
possible. The method we use is possibly the simplest one to implement. We remark that 
a better bound may be obtained by a more computationally intensive minimization. 

The method of reduction is simply a parametric reduction of the number of real fields 
in the Lagrangian of theory to 1. For instance one can define <pj = <x,$, j = 1, 2, n with 
Y^, otj — 1. This reduces a Lagrangian of n real fields to a Lagrangian of a single real field 
$. The advantage is an enormous simplification in the numerical search for the minimum 
of the action on C. Indeed, with spherical symmetry = Q(R),R 2 = x\ + ...x% the 
minimization is equivalent to solving the equation of motion: 

d 2 $ _ dU (d - 1) d$ 

dj¥~d$ R~dR ^'^ 

with the boundary conditions d^/dR = at R = and $ — > $^ as i? — > oc. There is a 
quick numerical method called the shooting method for solving this. One starts with a 
guess for the initial point $(0) and numerically integrates till either d^/dR changes sign 
(undershoot) or $ attains the value (overshoot). The initial point for the bounce must 
lie between the initial point of an overshoot and an undershoot, so the search converges 
rapidly by bisecting the interval between two initial point guesses of different outcomes. 

The following points in the method of reduction merit a special mention, 
(i) When some fields are put equal to zero (a, = 0) some coupling constants may drop out 
of the Lagrangian. The corresponding bounds on the parameter space are independent 
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of these couplings. Thus one can explore particular sets of parameters in the MSSM 
independently of the choice of other parameters. 

(ii) In general the potential U gives a good indication of the direction in the field space 
where the unstable vacuum is most likely to tunnel to. Therefore one can choose the 
reduction so that the reduced field $ points in the most obviously "dangerous" direction. 



3 The Unbounded From Below Directions in MSSM 

Tunneling to a generic CCB vacuum in the MSSM is difficult to explore using the for- 
malism of section 2 because of the presence of too many masses and couplings. We have 
chosen the directions in the field space where the potential may be unbounded from below 
(UFB) to compute the parameter bounds. In these directions, only a handful of the scalar 
fields are important and the reduced action contains fewer parameters. 

A general classification of UFB directions in the MSSM can be found in ref. |4[]. The 
unboundedness occurs along directions where the stabilizing D and F terms vanish and 
the soft terms are negative. This can happen in two different ways: 

(I) If Hi,H 2 and Lj are non-zero and all other scalar fields are zero. The index i is a 
generation index here. 

(II) If either 

(II a) Ei, Li, Lj and H 2 are non-zero and all other scalar fields are zero, or 
(II b) If Di, Qi, Lj and H 2 are non-zero and all other scalar fields are zero. 

The situations (II a) and (II b) lead to similar reduced potentials and will be treated 
together. Note that the effective potential is unbounded from below due to the tree 
approximation. The complete effective potential is bounded from below. However, when 
the quantum corrections to the tree order potential are small, tunneling amplitude to the 
deep CCB vacuum that replaces the UFB direction in the full effective potential should 
be calculable from the bounce computed using the tree order potential. This is because 
the bounce action actually depends on the "height" and "width" of the potential wall 
surrounding a false vacuum and is insensitive to the depth of the true vacuum. 

Let us write down the potential one obtains by setting all scalar fields equal to zero 
except the neutral Higgs fields H®, H 2 (where we write Hi = (H®, Hi), H 2 = (H 2 , H 2 )). 
A viable MSSM vacuum is obtained by minimizing this potential. 

U(H U H 2 ) = (ml i + \ l i\ 2 )Hl + (ml 2 +^)Hl-2mlH l H 2 + \[gl + gl] \H 2 2 - H 2 ] 2 (3.9) 
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where Hi = \Hf\,m^ = \Bfi\ and g± and g 2 are the U(l)y and SU{2) W gauge couplings 
respectively. We have chosen the phases of the fields to minimize the potential. Defining 
m 2 = m 2 H . + |/i| 2 , (i = 1, 2), we obtain the following parametric expressions for the vacuum 
expectation values (VEV's) v\ and v 2 for Hi and H 2 respectively. 

— = tan/5 (3.10) 

V2 



v. 



2 



2(motan 2 /5 — ml) 



2 ' g 2 {l -tan 2 /?) 

where g 2 = g\ + g\. This potential may have a local minimum at H\ = v 1, H 2 = v 2 with 
v 1 ,v 2 ^ if 

ml + ml > 2m 2 (3.11) 

2 2 4 

m 1 m 2 < m 3 . 

Either m 2 or m 2 can be negative. To get the correct W and Z masses we need 

v\ + v 2 w (175GeV) 2 (3.12) 

which fixes the value of to be 

2 _ |(mf + ml)|^/(m; + o)(m| + o) 
ms " |K + m 2 + 2a] I 

with a « ^ x (175GeV) 2 . 
3.1 Case I. 

Writing Lj = (z>^, the F terms are zero if H®, H® and Pi are non-zero but all other 
scalar fields are zero. The relevant scalar potential is: 



U(H, v) = m x 2 H\ + m\H\ - 2m\H x H 2 + m 2 v v 2 + ^g 2 [H 2 - H'f - v 2 \ ' (3.14) 

where v = \vl%\. Once again we have chosen the phases to minimize the potential. Notice 

that there are no F terms other than those that are already in equation ( |3.9| ). The D 

term vanishes identically along the surface H 2 — Hf — v 2 = 0. If we define y = H\/H 2 , 
then the potential on this surface is of the form 

U(H 2 ) = f(y)H 2 2 (3.15) 
f(y) = (ml + ml) - 2m\y + {m\ - ml)y 2 



2 
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The potential is unbounded from below if f(y) is negative. Minimizing f(y) we find the 



;; ' — {m{-ml) 



direction of fastest decreasing potential along y = y 

m\ — m 2 v > . 
To have a non-zero value of v 2 and H 2 we must have 

< y m < 1 • 

Finally, to have the unboundedness from below we require that 



provided 



fiVm) = (m 2 2 + m 2 u ) - 



mn 



'm 2 — m 2 ) 



< 0. 



(3.16) 



(3.17) 



(3.18) 




H 



( v l >*> ) 



Fig. 1 The reduced field. The piecewise straight line marked with arrows in the 
Hi,H 2 ,u plane defines a single real field. The line passes through the 
vacuum point H\ = v±, H2 = V2, v = 0. 



The last equation is the sufficiency constraint that we wish to relax to a necessary one. 
To compute the decay width of the viable vacuum we reduce the field space as follows. We 
take a piecewise straight line in the 3-dimensional space of Hi, H 2 and v parametrized by 
the real number e [—00, +00]. In the interval [—00, 0], we set |0| 2 = H\ + H\ with the 
points Hi,H 2 lying on the line connecting the point Hi = H 2 = v = and the (viable) 
vacuum point Hi = v i: H 2 = v 2 , v = 0. In the interval [0, 00], we define |0| 2 = Hf+H^+is 2 
with the points Hi,H 2 ,v lying on the line that passes through Hi = H 2 = v = and is 
defined by the relations Hi = y m H 2 , H\-H\-v 2 = 0, Hi, H 2 , v > [Fig. 1]. Along this 
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piecewise straight line the scalar potential reduces to 

f a^ + bcf 2 for <f) < 
U{(j>) = { (3.19) 

{ \f{y m ?<t> 2 for <p > . 

With this reduction the problem is reduced to solving a differential equation for a single 
real scalar. Note that the kinetic term for the field is T = / d 4 xJ2i(§^-) 2 which reduces 
to 27r 2 / RdR(^) 2 upon imposing 0(4) symmetry. We solve the differential equation for 
the bounce, on a computer, using the shooting method. Because the trilinear and Yukawa 
couplings do not enter any of the expressions in this case, the results obtained are valid 
for all three generations. 

3.2 Case II. 

Let us develop the formalism by considering the case (II a). In this case we have 
H 2 , E Ri , ELi,Vj 7^ 0, (i ^ j), and all other scalar fields are vanishing. We have retained 
the tilde-free notation of the previous case and set E^ = E Ri ,E Li = E Li . After setting 
F = \fiH 2 + h Ei E Li E Ri \ 2 = 0, and \E Ri \ = \E Li \ = E, one gets the potential: 

U(H 2 , E, v) = (mj - \^\ 2 )H 2 + (m 2 u + m 2 Et )E 2 + m%v 2 + \g 2 [h 2 + E 2 - v 2 \ , (3.20) 

where we have simplified notation with \vj\ —>■ v. Note that h E i is just a number. The D 
term is zero if v 2 = H\ + E 2 . The potential in this direction is: 

U'(H 2 ) = (mj - |/i| 2 + m 2 Lj )H 2 ± M( m ^ + m |. + m %)H 2 . (3.21) 

n>Ei 

The ambiguous sign comes from phase choices of the fields. We will choose the sign that 
makes the second term on the right negative (i.e. a faster drop in the potential). The 
resulting sufficiency condition for vacuum stabilty is a strong constraint (the Komatsu 
constraint 0): 

m 2 2 - \fi\ 2 + m 2 Lj > 0. (3.22) 

The field reduction is now slightly more complicated than that of the previous case. We 
define a real scalar G [— oo,+oo] that for cf) > is the path length of the curved 
line defined by the equations H\ + E 2 — v 2 = and F = measured from the point 
H 2 = v = E = 0. When (j) < it is given by <\> 2 = Hf + R\ with Hi = y m H 2 , as in 
the previous section. The infinite (curved) line <fi G [— oo, +oo] lies in the 5- dimensional 
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Euclidean space with coordinates Hi, H 2 , Eu, Em, v. The scalar potential along this 
curved line is 

f a0 4 + b<p 2 for < 
U(<f>) = (3.23) 
{ U'{H 2 ) for 0>O. 

The parametric representation of the potential in terms of H 2 hides the main complexity of 
this case, namely, (when positive) must be computed as a function of H 2 by integrating 
a differential along the curve defined above: 



[dH$ + dE 2 Li + dE 2 Ri + dv 2 } . 



(3.24) 



The case II b is trivially obtained by the substitution E { — > Di, Li — > Qi and h Ei — > h Di 
in equations ( p.20| ) and (|3.21|) . 

4 Results 
4.1 Case I. 




1 10 100 

M 4 

2 



Fig. 2. Bounce actions in Case I. We have drawn iso-action lines connecting 
points with action (shown adjacent to lines) of same order. All masses 
are in units of 100 GeV. The regions marked "I" on either side of broken 
lines have no points relevant to tunneling computation. 



In this case the requirement that a viable MSSM vacuum exist (equations ( |3.11| , |3.13|) ) 
is hard to satisfy simultaneously with the requirement that destabilizing directions exist 
(equations ( |3.16| , |3.17] , |3.18|) ). It is convenient to plot the results in the Mf,M% plane 
where we define Mf = (m 2 + m 2 u )(m\ — m 2 ) and M 2 = m\. The compression of the 4- 
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dimensional parameter space (mi, m 2 , /i, mu) to this plane leads to little loss of infor- 
mation since the bounce action is relatively insensitive to the other two axes. In Fig. 2 
regions marked with an "I" do not satisfy the requirements mentioned above. This leaves 
a narrow window in the realistic mass range open for our computations where the bounce 
actions are greater than 400h (usually by an order of magnitude). We have drawn iso- 
action lines for the bounces to display the action variation. Since we are unable to rule out 
any part of the parameter space with certainty, this case serves to display the limitation 
of our algorithm. However the plot shows the range and variation of the bounce action 
and the result here may be improved by a more accurate numerical calculation. 



4.2 Case II a. 

We have reduced the number of free parameters by using approximate degeneracy of the 
slepton masses. The restriction to this subspace of the parameter space is motivated by 
Grand Unification of couplings and masses (for a review see ref. 0). To be conservative, 
we chose the case of the r slepton (i.e. i = 3 in the relevant formulas) which has the 
largest Yukawa coupling and (given slepton mass unification) from inspection of equation 
( p.21| ), is likely to give the least negative contribution to U'(H 2 ). The results are not 
very sensitive to slepton mass splittings. Again it is convenient to define the masses 
M| = — (m 2 , — |/i| 2 + an d — ^~\^\ m %- The plot of bounce actions on the 

M| vs. Mf plane is relatively insensitive to other independent axes and shows a large 
region where our methods are adequate to prove the instability of the MSSM vacuum 
(Fig. 3). We have scanned more than 2 x 10 5 points in the the parameter region m 2 G 
[-100, 500], m\ G [-2, 100], /i 2 G [0,500] and m 2 Lj G [0.2,100] in units of (100 GeV) 2 . A 
fraction of these points appear in Fig. 3 which contains our result in the form of the broken 
line. To the right of the line a stable and viable MSSM vacuum is ruled out because at 
least one of the following necessarily happens: (i) equations ( |3.11|) , (|3.13j ) are not satisfied, 



(ii) slepton mass < 45 GeV, (iii) | < 400. To the left of the line a stable and viable MSSM 
vacuum may or may not exist. To obtain this line, one need only plot those points for 
which equations ( |3.11| ,), ( |3.13| ) are satisfied and slepton mass < 45 GeV and then draw 



the line to bound all points with ^ > 400 (the crosses) from the right. Indeed, the points 
with f < 400 need not be plotted either. In Fig. 3, we plot a fraction of the points that 
have ^ < 400 to illustrate our coverage of the parameter space. 

The broken line is therefore our main result. It is a modification of the Komatsu 
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constraint to the approximate relation for vacuum instability: 



Mix + 0.2MX > 9.5x d 



(4.25) 



valid for M 3 > with x = 100 GeV and when all scalar masses in the theory are within two 
orders of magnitude of the electroweak scale. The simple method used by us is successful 
in ruling out most of the relevant region in the parameter space. 
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Fig. 3. Bounce actions in Case II a. Crosses denote points with S/% > 400. The 
circles mark points with S/fi < 400. The action tends to decrease as one 
moves to the right or upwards of the broken line and only points with 
S/% < 400 appear outside the window we display. All masses are in units 
of 100 GeV. 



4.3 Case II b. 

This case is identical to the previous one with the replacement L; L — > Q^E i — * Di and 
hEi — ► hoi- Again we use approximate degeneracy of squark and slepton masses and 
study the bottom squark which has the largest Yukawa coupling (i = 3). We define 
M| = — (ml — |/i| 2 +iriQj), and Mf = j^\fi\m 2 Lj . Scanning the same parameter region as 
in the previous case we obtain a plot of bounce actions on the Mf vs. Mf plane (Fig. 4) 
which is similar to the plot in the previous case and rules out a comparable volume in 
the parameter space. The corresponding modification of the Komatsu constraint is the 
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following relation for vacuum instability: 



Mix + 0.15M? > 6.0a; 



(4.26) 



valid for M 5 > 0, with x = 100 GeV. 
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Fig. 4. Bounce actions in Case II b. The conventions of Fig. 3 apply. 



As mentioned in section 3, the action values are rigorous upper bounds on the bounce 
action provided one rescales the bounce points to lie on the space C. The rescaling factor 
departs from 1 by less than 5 percent ([T + 2V] 2 /T 2 is less than 10 -3 in all cases). We 
also found the bounce size to be between 1 and 10 in units of (100 GeV) -1 which justifies 
the approximation rj — 1. 



5 Conclusions 

In conclusion, we have used a new and simple method of computing vacuum decay rates 
in field theories to scan the parameter space of the MSSM and test for vacuum stability. 
We find that in one class of the characteristic UFB direction, our method fails to yield a 
useful constraint. However it does successfully compute a lower bound on the decay rate. 
In another class of the UFB direction our method yields very strong constraints and we 
are successful in ruling out with certainty a significant part of the parameter space where 
the MSSM vacuum is met ast able. 
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